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Abstract The MERRA Aerosol Reanalysis (MERRAero) has been recently developed at NASA's Global 
Modeling Assimilation Office. This reanalysis is based on a version of the Goddard Earth Observing System-5 
(GEOS-5) model radiatively coupled with Goddard Chemistry, Aerosol, Radiation, and Transport aerosols, and it 
includes assimilation of bias-corrected aerosol optical thickness (AOT) from the Moderate Resolution Imaging 
Spectroradiometer (MODIS) sensor on both Terra and Aqua satellites. In October over the period 2002-2009, 
MERRAero showed that AOT was lower over the east of the Ganges basin than over the northwest of the 
Ganges basin: this was despite the fact that the east of the Ganges basin should have produced higher 
anthropogenic aerosol emissions because of higher population density, increased industrial output, and 
transportation. This is evidence that higher aerosol emissions do not always correspond to higher AOT over 
the areas where the effects of meteorological factors on AOT dominate those of aerosol emissions. MODIS 
AOT assimilation was essential for correcting modeled AOT mainly over the northwest of the Ganges basin, 
where AOT increments were maximal. Over the east of the Ganges basin and northwest Bay of Bengal 
(BoB), AOT increments were low and MODIS AOT assimilation did not contribute significantly to modeled 
AOT. Our analysis showed that increasing AOT trends over northwest BoB (exceeding those over the east of 
the Ganges basin) were reproduced by GEOS-5, not because of MODIS AOT assimilation but mainly because of the 
model capability of reproducing meteorological factors contributing to AOT trends. Moreover, vertically integrated 
aerosol mass flux was sensitive to wind convergence causing aerosol accumulation over northwest BoB. 


1. Introduction 

The Indian subcontinent (and the Ganges basin in particular) is characterized by a significant population 
growth accompanied by developing industry, agriculture, and increasing transportation. This has resulted in 
declining air quality [Di Girolamo et a!., 2004; Ramanathan and Ramana, 2005; Tripathi et al., 2005; Prasad and 
Singh, 2007; Kaskaoutis et al., 2011a; Dey and Di Girolamo, 201 1; Krishna Moorthy et al., 2013]. With respect to 
air pollution, one could suggest some relationship between population figures and anthropogenic aerosol 
emissions. Indeed, Kishcha et al. [2011] showed that over extensive areas with differing population densities 
in the Indian subcontinent, the higher the averaged population density, the larger the averaged aerosol 
optical thickness (AOT). In addition, the larger the population growth, the stronger the increasing AOT trends. 

In accordance with Di Girolamo etal. [2004], Prasad and Singh [2007], and Kumar etai. [2010], prevailing winds 
blowing along the Ganges basin in the postmonsoon and winter months transport anthropogenic aerosol 
particles into the Bay of Bengal. A number of sea expeditions to BoB were conducted to investigate the 
resulting increased levels of air pollution over BoB [Ramachandran and Jayaraman, 2003; Vinoj et al., 2004; 
Ganguly et al., 2005; Moorthy et al., 2008; Kumar et al., 2010; Kaskaoutis et al., 201 1 b]. Moreover, long-term AOT 
trends over South Asia, including BoB, were examined, using different satellite AOT data sets, by Mishchenko 
and Geogdzhayev [2007], Zhao et al. [2008], Zhang and Reid [201 0], Kaskaoutis et al. [201 1 a], Dey and Di 
Girolamo [2011], and Hsu et al. [2012]. Based on advanced very high resolution radiometer (AVHRR) satellite 
data, Mishchenko and Geogdzhayev [2007] compared overwater AOT averaged over two separate periods, 

1 988-1 991 and 2002-2005, and found significant changes. Zhao et al. [2008] studied AOT trends over the 
whole area of BoB for spring, summer, autumn, and winter during the 25 year period 1 981-2005, using 
AVHRR data. Using Moderate Resolution Imaging Spectroradiometer (MODIS)-Terra level 2 AOT data, Zhang 
and Reid [2010] analyzed AOT trends over the whole area of BoB for all months during the 10 year period 
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2000-2009. The spatial distribution of decadal (2000-2009) MODIS level-3 AOT trends over South Asia, 
including BoB, in different months was obtained by Kaskaoutis et al. [201 la]. Using Multiangle Imaging 
SpectroRadiometer aerosol data, decadal (2000-2009) AOT trends over the Indian subcontinent and surrounding 
sea areas were also estimated by Dey and Di Girolamo [201 1]. Hsu etal. [2012] created maps of Sea-viewing Wide 
Field-of-view Sensor AOT trends over the period 1998-2010 for each of the four seasons. 

The early postmonsoon season over the study region is characterized by aerosol transport from the Ganges 
basin to northwest BoB by prevailing winds and still significant rainfall of over 150 mm/month over the east 
of the Ganges basin and northwest BoB. It would be reasonable to consider that AOT trends over sea areas in 
BoB were created by changes in aerosol sources on the land in the Indian subcontinent. In our previous study 
[Kishcha et a!., 2012], we found that it was not always the case. Specifically, we found that in October, MODIS 
showed strong increasing aerosol optical thickness (AOT) trends over northwest Bay of Bengal (BoB) in the 
absence of AOT trends over the east of the Indian subcontinent. This was unexpected, because sources of 
anthropogenic pollution were located over the Indian subcontinent, mainly in the Ganges basin, and aerosol 
transport from the Indian subcontinent to northwest BoB was carried out by prevailing winds. 

It was our purpose to determine whether existing state-of-the-art aerosol data-assimilated systems were 
capable of reproducing the aforementioned AOT trends over northwest BoB in the early postmonsoon season, 
in the presence of significant rainfall. For the model, it would be a challenge just to obtain correct space-time 
distribution of rainfall, which is of importance for estimating aerosol wet removal. The NASA Goddard Earth 
Observing System-5 (GEOS-5) was used to extend the NASA Modern Era-Retrospective Analysis for Research 
and Applications (MERRA) reanalysis by adding five atmospheric aerosol components (sulfates, organic carbon, 
black carbon, desert dust, and sea salt). In the current study, the obtained 8 year (2002-2009) assimilated 
aerosol data set (so-called MERRAero) was applied to examine aerosol trends over the Ganges basin and 
northwest Bay of Bengal (BoB) in the postmonsoon season. In addition, using air pollution modeling allowed 
us to look at the situation over the study region from the point of view of air quality. This was carried out by 
estimating the contribution of various aerosol species to AOT and its trends. 

Note that AOT assimilation was effective only for two short periods of MODIS's appearance over the study 
area. All other times (1 8 h per day) the GEOS-5 model worked independently of MODIS. 

2. GEOS-5 and the MERRA Aerosol Reanalysis (MERRAero) 

2.1. GEOS-5 Earth Modeling System 

GEOS-5 is the latest version of the NASA Global Modeling and Assimilation Office (GMAO) Earth system model. 
GEOS-5 contains components for atmospheric circulation and composition (including atmospheric data 
assimilation), ocean circulation and biogeochemistry, and land surface processes. Components and individual 
parameterizations within components are coupled under the Earth System Modeling Framework [Hill 
et al., 2004], In addition to traditional meteorological parameters (winds, temperatures, etc.) [Rienecker 
et al., 2008], GEOS-5 includes modules representing the atmospheric composition, most notably aerosols 
[Colarco et al., 2010], and tropospheric/stratospheric chemical constituents [ Pawson et al., 2008], and the 
impact of these constituents on the radiative processes of the atmosphere. 

2.2. Aerosols in GEOS-5 

GEOS-5 includes modules representing atmospheric composition, including aerosols [ Colarco etal., 2010] and 
tropospheric and stratospheric chemical constituents [Pawson et al., 2008]. The current generation aerosol 
module is based on a version of the Goddard Chemistry, Aerosol, Radiation, and Transport (GOCART) model 
[Chin et al., 2002]. GOCART treats the sources, sinks, and chemistry of dust, sulfate, sea salt, and black and 
organic carbon aerosols. Aerosol species are assumed to be external mixtures. Both dust and sea salt have wind 
speed-dependent emission functions [Colarco et al., 2010], while sulfate and carbonaceous species have 
emissions principally from fossil fuel combustion, biomass burning, and biofuel consumption, with additional 
biogenic sources of organic carbon. Sulfate has additional chemical production from oxidation of S0 2 
and dimethylsulfide and a database of volcanic S02 emissions and injection heights. Aerosol emissions 
for sulfate and carbonaceous species are based on the AeroCom version 2 hindcast inventories (T. Diehl, 
personal communication, 2013) (http://aerocom.met.no/emissions.html). Daily biomass burning emissions 
are from the Quick Fire Emission Dataset (QFED) and are derived from MODIS fire radiative power retrievals 
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Table 1. The Overview of Main Attributes of NASA MERRAero Assimilated Aerosol Data 


Feature 


Description 


Model 


GEOS-5 Earth Modeling System (with GOCART aerosol components) 
Constrained by MERRA Meteorology (Replay) 

Land sees observed precipitation (like MERRA Land) 

Driven by QFED daily Biomass Emissions 
Local displacement ensembles 
MODIS reflectances 

AERONET 3 calibrated AOT's (neural net) 

Stringent cloud screening 
Middle 2002 to present (Aqua + Terra) 

Florizontal: nominally 50 km 
Vertical: 72 layers, top ~85 km 
Dust, sea salt, sulfates, and organic and black carbons 


Aerosol data assimilation 


Period 

Resolution 


Aerosol species 


a AERONET, Aerosol Robotic Network. 

[Darmenov and da Silva, 2013]. Total mass of sulfate and carbonaceous aerosols are tracked, while for dust 
and sea salt the particle size distribution is explicitly resolved across five noninteracting size bins for each. 

For all aerosol species, optical properties are primarily from the commonly used Optical Properties of Aerosols 
and Clouds (OPAC) data set [Hess et a/., 1998]. OPAC provides the spectrally varying refractive index and a 
humidification factor for each aerosol species which, together with assumptions about the particle size dis- 
tribution of each species, are used to construct spectrally varying lookup tables of aerosol optical properties 
such as the mass extinction efficiency, single scattering albedo, and asymmetry parameter, inputs required 
by our radiative transfer codes (details are in Colarco et at. [2010, and references therein]). 

2.3. GEOS-5 Data Assimilation 

GEOS-5 has a mature atmospheric data assimilation system that builds upon the Gridpoint Statistical 
Interpolation (GSI) algorithm, jointly developed with National Centers for Environmental Prediction (NCEP) 
[Wu et at., 2002; Derber et at., 2003; Rienecker et al., 2008]. The GSI solver was originally developed at NCEP as 
an unified three-dimensional variational analysis system for supporting global and regional models. GSI in- 
cludes all the in situ and remotely sensed data used for operational weather prediction at NCEP. 

GEOS-5 also includes assimilation of AOT observations from the MODIS sensor on both Terra and Aqua sat- 
ellites. Based on the work of Zhang and Reid [2006] and Lary et al. [2009], a back-propagation neural network 
has been developed to correct observational biases related to cloud contamination, surface parameteriza- 
tion, aerosol microphysics, etc. Online quality control is performed with the adaptive buddy check of Dee etal. 
[2001], with observation and background errors estimated using the maximum likelihood approach of Dee 
and da Silva [1999]. 

2.4. MERRA Aerosol (MERRAero) Reanalysis 

MERRA is a NASA reanalysis for the satellite era using a major new version of the Goddard Earth Observing 
System Data Assimilation System Version 5 (GEOS-5) [Rienecker et al., 201 1], The project focuses on historical 
analyses of the hydrological cycle from the NASA EOS suite of observations in a climate context, on a broad 
range of weather and climate time scales and places. The MERRA time period covers the modern era of re- 
motely sensed data, from 1 979 to the present, and the special focus of the atmospheric assimilation is the 
hydrological cycle. Like other similar reanalysis, MERRA provides meteorological parameters (winds, tem- 
perature, and humidity), along with a number of other diagnostics such as surface and top of the atmosphere 
fluxes, diabatic terms, and the observational corrections imposed by the data assimilation procedure. 

As a step toward an Integrated Earth System Analysis, the GMAO is producing several parallel reanalyses of 
other components of the Earth system such as ocean, land, and atmospheric composition. Of particular rel- 
evance for this paper, the MERRA Aerosol Reanalysis (MERRAero), where MODIS AOT observations are as- 
similated providing a companion aerosol gridded data sets that can be used to study the impact of aerosols 
on the atmospheric circulation and on air quality in general. Table 1 summarizes the main attributes of 
MERRAero. Notice that MERRAero only covers the later years of MERRA, capitalizing on the improved aerosol 
measurements from NASA's EOS platforms. 
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Figure 1. Spatial distributions of (a) the 8 year (2002-2009) mean MERRAero AOT and (b) its trends (characterized by AOT slopes) in 
October. The AOT trend values correspond to the slope of the linear regression analysis. The squares show the locations of zones 1 to 
7 within the study region. 


3. Method 

Following our previous study [Kishcha eta!., 2012], we analyzed long-term variations of AOT over seven zones, 
each 3° x 3°, located in the Ganges basin and northwest BoB (Figure 1 ). As mentioned, in the postmonsoon 
period, prevailing winds blow along the Ganges basin. The specified zones in the Ganges basin provide us 
with an opportunity for analyzing air pollution trends produced by local sources and aerosol transport. 
Figure la shows the spatial distribution of 8 year mean MERRAero AOT over the region under consid- 
eration in October, together with the location of zones 3° x 3° in the Indian subcontinent (zone 1 to 
zone 5) and in the Bay of Bengal (zones 6 and 7). MERRAero monthly AOT data are available from the year 
2002. To analyze AOT and its trends over the Indian subcontinent and northwest BoB, we used monthly 
MERRAero AOT data with horizontal resolution of approximately 50 km, during the 8 year period 2002-2009. To 
estimate the NASA GEOS-5 model performance, 3 h MERRAero AOT data together with AOT increments 
were used. 

A linear fit was used to determine the resulting trend of aerosol optical thickness during the study period 
(2002-2009) over each of the aforementioned zones. The obtained AOT trend values correspond to the slope 
of the linear fit. To ensure that the linear fit produced normally distributed residuals, they were required to 
pass the Shapiro-Wilk normality test [Shapiro and Wilk, 1965; Razali and Wah, 2011]. If the residuals were 
normally distributed, they could be used in a f test, in order to estimate the statistical significance of a linear 
fit. The statistical significance of the AOT trend was checked by applying the significance level (p) value, i.e., 
p < 0.05 for statistically significant AOT trends at the 95% confidence level. 

4. Results 

4.1. Total MERRAero AOT and Its Trends in October 

In accordance with spatial distribution of 8 year mean AOT in the early postmonsoon season (October), 
MERRAero showed high AOT values over the Ganges basin with a maximum over the northwest part of the 
Ganges basin (Figure la). Therefore, MERRAero data were able to reproduce the main structure of aerosol 
distribution over the Ganges basin. The Ganges basin is the most polluted part of the Indian subcontinent, 
where highly populated areas and main industrial centers are located. 

We analyzed zone-to-zone variations of MERRAero AOT averaged over the specified zones. In the early 
postmonsoon season (October), MERRAero showed mainly decreasing AOT variations from zone 3 to zone 5 
(Figure 2a and Table 2). Note that this decrease in AOT from northwest to east of the Ganges basin does not 
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correspond to the distribution of population 
density: population density is higher in the east of 
the Ganges basin (zones 4 and 5) than in the 
northwest of the Ganges basin (zone 1) (Figure 3). 
At first glance, this is contradictory to our previous 
findings on the relationship between AOT and 
population density in the Indian subcontinent 
[Kishcha et at., 201 1], It should be mentioned, 
however, that in our previous study, we used 
averaging over significant areas of the Indian 
subcontinent with differing population densities. 

zone 1 zone 2 zone 3 zone 4 zone 5 zone 6 zone 7 



The most probable reason for the decrease in AOT 
over the east of the Ganges basin, where popu- 
lation density is the highest in the entire Ganges 
basin, is wet removal processes after significant 
rainfall in October. Monthly accumulated Tropical 
Rainfall Measuring Mission (TRMM) rainfall data 
from the 3B42V6 archive, on a 0.25° x 0.25° lati- 
tude-longitude grid [Huffman et at., 2007], were 
used to estimate zone-to-zone variations of 8 year 
(2002-2009) mean TRMM rainfall over the speci- 
fied zones in October (Figure 4a). High rainfall 
values of over 150 mm can be seen in October 
over the east of the Ganges basin (zone 5) and 
northwest BoB (zone 6). Moreover, rainfall data 
showed that over the east of the Ganges basin, 
the accumulated rainfall in October in the first 4 
year period 2002-2005 was essentially higher 
than in the second 4 year period 2006-2009 

(Figure 4a). As a result, higher values of MERRAero AOT over the east of the Ganges basin (zones 4 and 5) were 
observed in the second 4 year period 2006-2009 than in the first 4 year period 2002-2005 (Figure 2b). 
Moreover, the aforementioned decrease in the 8 year (2002-2009) mean AOT over the east of the Ganges basin 
and northwest BoB in October corresponds to the increase in model-simulated aerosol wet deposition over the 
specified zones toward northwest BoB, in accordance with the spatial distribution shown in Figure 4b. 



Figure 2. (a) Zone-to-zone variations of 8 year (2002-2009) mean 
MERRAero AOT averaged over the specified zones, (b) Zone-to-zone 
variations of MERRAero AOT averaged over the first 4 year (2002-2005) 
period and over the second 4 year (2006-2009) period. The error bars 
show the standard error of mean AOT. 


Therefore, we have arrived at the interesting point, which is as follows. The east of the Ganges basin should 
have produced higher anthropogenic aerosol emissions than the northwest of the Ganges basin because of 


Table 2. Eight Year (2002-2009) Mean AOT (r), Standard Deviation (SD), and AOT Slope (a) of MERRAero AOT Averaged Over 
the Specified Zones in October 3 


Area 

Zone # 

Geographic Coordinates 

r 

SD 

a (Per Year) 

S-W Test 

P 

IS 

1 

28.5°N-31.5°N 

72.7°E-75.7°E 

0.53 

0.07 

-0.005 

Normal 

Not significant 


2 

27°N-30°N 

75.7°E-78.7°E 

0.50 

0.08 

-0.004 

Normal 

Not significant 


3 

25°N-28°N 

78.7°E-81.7°E 

0.51 

0.08 

0.008 

Normal 

Not significant 


4 

24°N-27°N 

82.5°E-85.5°E 

0.43 

0.05 

0.012 

Normal 

Not significant 


5 

22°N-25°N 

86°E-89°E 

0.35 

0.05 

0.010 

Normal 

Not significant 

BoB 

6 

18°N-21°N 

87°E-90°E 

0.22 

0.05 

0.015 

Normal 

0.043 


7 

15°N-18°N 

84°E-87°E 

0.20 

0.05 

0.020 

Normal 

0.006 


a The decision based on the Shapiro-Wilk normality test for residuals (S-W test) and the significance level (p) are also displayed. If the 
p value was too high as compared with the 0.05 significance level, the obtained linear fit was considered as statistically insignificant. 
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higher population density, 
increased industrial output and 
transportation. Nevertheless, 
MERRAero showed that AOT over 
the east of the Ganges basin was 
lower than that over the northwest 
of the Ganges basin. This is 
evidence that higher aerosol emis- 
sions do not always correspond to 
higher AOT. This can take place 
over the areas where the effects of 
meteorological factors on AOT (in 
this case, precipitation) dominate 
those of aerosol emissions. 

Spatial distributions of MERRAero 
AOT trends during the 8 year (2002- 
2009) study period showed strong 
increasing AOT trends over north- 
west BoB exceeding those over the 
Ganges basin (Figure 1b). This indi- 
cates that MERRAero is capable of 
reproducing the main features of 
the phenomenon of strong increas- 
ing AOT trends over northwest BoB 
in the early postmonsoon season, in 
line with our previous study [Kishcha 
et at., 2012], 

4.2. Effects of Rainfall on 
MERRAero AOT 

As mentioned, in the early 
postmonsoon season (October), 
intense rainfall can be frequently 
observed over the east of the 
Ganges basin. These severe pre- 

Figure 3. (a) Population density (persons km ) distributions over the Indian subcon- 
tinent. (b) Zone-to-zone variations of population density averaged over the specified Cipitation events COLlId Strongly 
zones. The GPW-v3 gridded population density data for the year 2005 were used affect AOT Over the east of the 
(http://sedac.ciesin.coiumbia.edu/data/coiiection/gpw-v3). Ganges basin due to aerosol wet 

removal processes. To understand 

the rain effects on AOT over the east of the Ganges basin (zone 5), we compared year-to-year variations of 
assimilated MERRAero AOT and TRMM accumulated rainfall, over zone 5 in each October during the study 
period (not shown). Rainfall data showed that the accumulated rainfall in October in the first 4 year period 
2002-2005 was higher than in the second 4 year period 2006-2009. A strong inverse relationship (with a high 
negative correlation of over —0.8) between changes in assimilated MERRAero AOT and rainfall is clearly seen: 
each increase in rainfall was accompanied by a decrease in assimilated AOT. The aforementioned decrease in 
rainfall over zone 5 in October during the study period can explain some increasing trend in MERRAero AOT 
observed over that area in October (Table 2). There was some dissimilarity in the rainfall amount between the 
east of the Ganges basin (zone 5) and northwest BoB (zone 6): northwest BoB does not show as clear decreasing 
trends in rainfall amount as the east of the Ganges basin does. 

4.3. AOT of Different Aerosol Species and Their Trends in the Early Postmonsoon Season 

As known, satellite remote sensing data cannot distinguish between various aerosol species. MERRAero pro- 
vides us with an opportunity to look at the situation over the study region from the point of view of air quality. 
This was carried out by estimating the contribution of various aerosol species to AOT and its trends. Based on 
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Figure 4. (a) Zone-to-zone variations of TRMM accumulated rainfall over the specified 
zones in October averaged over the 8 year study period (2002-2009), over the first 4 year 
period (2002-2005), and over the second 4 year period (2006-2009). The error bars show 
the standard error of mean accumulated rainfall. TRMM data from the 3B42V6 archive were 
used, (b) The spatial distribution of the 8 year (2002-2009) mean model-simulated aerosol 
wet deposition rate (kg m 2 s ') in October. 


MERRAero model data, Figure 5a 
represents zone-to-zone variations 
of 8 year (2002-2009) mean AOT 
of several aerosol components 
(desert dust, organic and black 
carbon, and sulfates) averaged 
over specified zones in October, 
and their trends. One can see that 
over the west of the Ganges basin 
(zones 1-3), where precipitation 
was minimal (Figure 4a), each 
aerosol component changes 
from zone-to-zone in a different 
way, in accordance with its own 
emissions, transport, and gravi- 
tational settling (Figure 5a). By 
contrast, over the east of the 
Ganges basin (zones 4 and 5), a 
general decrease in AOT can be 
observed due to increasing 
aerosol wet removal toward 
northwest BoB (Figure 4b). Over 
zone 1, there is a considerable 
amount of carbon aerosols (as 
a result of crop waste burning 
[Sharma et at., 2010]), dust par- 
ticles, and sulfate aerosols 
(Figure 5a). This explains the 
AOT maximum over the north- 
west of the Ganges basin in 
October. Over sea areas (zones 
6-7), aerosols are dominated by 
anthropogenic air pollution, such 
as sulfates and carbon aerosols 
(Figure 5a). This predominance 
can affect cloud formation, 
atmospheric dynamics, and 
even marine life in this region. 


By contrast to sulfates and car- 
bonates, dust aerosol particles have no sources along the Ganges basin. Therefore, dust distribution 
along the Ganges basin is determined by aerosol transport (by the action of prevailing winds blowing 
along the Ganges basin) and by deposition processes. One can see that the 8 year mean dust AOT 
values noticeably decreased along the Ganges basin and over northwest BoB. This resulted in the 
decrease in dust contribution to the total AOT from approximately 30% over zone 1 to 8% over zones 
from 5 to 7 (Table 3). Furthermore, dust AOT trends did not change in transition from land to sea: ap- 
proximately the same slightly increasing dust AOT trends of -0.004 yr -1 were obtained along the east of 
the Ganges basin and over northwest BoB (Figures 5b and 6c). We found that these AOT trends over 
zones from 5 to 7 were statistically significant (Table 3). The same dust AOT trends along the Ganges 
basin and over northwest BoB suggest an increasing trend in some external source of dust emissions, 
outside the Ganges basin. It should be kept in mind that MERRAero only assimilates total AOT and that 
the trend in aerosol speciation may depend on the trend (or lack thereof) of the specified emissions. 


The distribution of sulfate AOT along the Ganges basin is determined by sulfate aerosol emissions, 
together with aerosol transport (by the action of prevailing winds) and deposition processes (Figure 5a). The 
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Figure 5. Zone-to-zone variations of (a) 8 year (2002-2009) mean MERRA AOT 
of various aerosol species (sulfates (SU), organic and black carbon (OC and BC), 
and desert dust) averaged over the specified zones in October, and (b) their 
AOT trends (characterized by the slope of the linear regression analysis). The 
error bars show the standard error of mean AOT. 


sulfate contribution to the total AOT increased 
along the Ganges basin from approximately 
30% over zone 1 to -56% over zone 5 (Table 3). 
Over northwest BoB (zones 6 and 7), the sulfate 
contribution to the total AOT was over 50% 
(Table 3). Thus, according to MERRAero AOT 
data, sulfates were the major atmospheric 
aerosol component over the east of the 
Ganges basin and over northwest BoB. 
Moreover, MERRAero data showed that sulfate 
AOT trends changed in transition from land to 
sea: strong statistically significant increasing 
sulfate AOT trends (of 0.008 yr _1 and 
0.01 1 yr~' over zones 6 and 7, respectively) 
exceeded those over the east of the Ganges 
basin (zone 5) (Figures 5b and 6a and Table 3). 

With respect to organic and black carbon 
aerosols, their distribution of 8 year mean 
AOT values along the Ganges showed a wide 
maximum from the northwest to the center 
of the Ganges basin (zones from 1 to 3) 
(Figure 5a). As mentioned, this area of maxi- 
mum carbon AOT is known for crop waste 
burning aerosols [ Sharma et at., 2010; 
Venkataraman et at., 2006]. AOT values of 
carbon aerosols decrease to the east from 
zone 3 (Figure 5a). As discussed in section 4.1, 
the reason for the decrease in AOT over the 


Table 3. The Eight Year Mean AOT (r), Standard Deviation (SD), and AOT Slope (a) for Long-Term Changes of MERRAero AOT 
for Different Aerosol Species (Desert Dust, Organic and Black Carbon, and Sulfates) Averaged Over Specified Zones in October 3 


Area 

Zone # 

F (%) 

r 

SD 

a (Per Year) 

S-W Test 

P 

IS 

1 

31.7 

0.17 

Sulfates 

0.04 

0.001 

Normal 

Not significant 


2 

36.7 

0.18 

0.04 

-0.001 

Normal 

Not significant 


3 

43.7 

0.22 

0.05 

0.004 

Normal 

Not significant 


4 

49.7 

0.22 

0.04 

0.006 

Normal 

Not significant 


5 

55.9 

0.20 

0.03 

0.004 

Normal 

Not significant 

BoB 

6 

52.9 

0.12 

0.03 

0.008 

Normal 

0.050 


7 

50.7 

0.10 

0.03 

0.011 

Normal 

0.004 

IS 

1 

37.7 

0.20 

Organic and Black Carbon 

0.05 -0.01 1 

Normal 

Not significant 


2 

39.2 

0.20 

0.03 

-0.009 

Normal 

Not significant 


3 

37.7 

0.19 

0.02 

-0.001 

Normal 

Not significant 


4 

37.5 

0.16 

0.02 

0.002 

Normal 

Not significant 


5 

34.9 

0.12 

0.02 

0.003 

Normal 

Not significant 

BoB 

6 

28.8 

0.06 

0.02 

0.004 

Normal 

Not significant 


7 

25.9 

0.05 

0.02 

0.005 

Normal 

0.026 

IS 

1 

29.5 

0.16 

Desert Dust 
0.04 

0.006 

Normal 

Not significant 


2 

23.2 

0.12 

0.03 

0.006 

Normal 

Not significant 


3 

17.8 

0.09 

0.03 

0.005 

Normal 

Not significant 


4 

12.2 

0.05 

0.02 

0.004 

Normal 

Not significant 


5 

8.0 

0.03 

0.01 

0.004 

Normal 

0.035 

BoB 

6 

8.3 

0.02 

0.01 

0.004 

Normal 

0.012 


7 

8.1 

0.02 

0.01 

0.004 

Normal 

0.008 


3 F corresponds to the fraction of aerosol component AOT (in percentages) from the total MERRAero AOT. 
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east of the Ganges basin in October is significant 
rainfall accompanied by aerosol wet removal 
processes. The joint contribution of organic and 
black carbon aerosols to the total AOT is -38% 
over the northwest of the Ganges basin (zones 
from 1 to 3), -35% over the east of the Ganges 
basin (zone 5), and approximately 27% over 
northwest BoB (zones 6 and 7) (Table 3). Similar to 
AOT trends of sulfate aerosols, MERRAero showed 
that AOT trends of carbon aerosols changed in 
transition from land to sea: increasing AOT trends 
in organic and black carbon AOT over the sea 
(zones 6 and 7) exceeded those over zone 5 in the 
land (Figures 5b and 6b and Table 3). 

Based on MERRAero data, we found that in 
October the contribution of sea-salt aerosols to 
the total AOT over the east of the Ganges basin 
was even lower than that of desert dust. Over 
northwest BoB, desert dust and sea-salt aerosols 
equally contributed to the total AOT (Figure 7a). 
Our analysis showed that sea-salt aerosols did not 
contribute to the increasing AOT trends over 
northwest BoB: no sea-salt AOT trend was ob- 
served in year-to-year variations during the study 
period (Figure 7a). This was supported by year-to- 
year variations of MERRA surface winds over 
northwest BoB (Figure 7b). 

4.4. Factors Contributing to AOT Trends Over 
Northwest BoB 

MERRAero showed increasing AOT trends over 
northwest BoB in October exceeding AOT trends 
over the east of the Ganges basin (Figure 1b). This 
was despite the fact that sources of air pollution 
are located on the land, mainly in the Ganges 
basin. There could be several factors contributing 
to the increasing AOT trends over northwest BoB. 
First, there were changes in the atmospheric cir- 
culation over northwest BoB in October during 
the 8 year study period (Figure 8). Mean wind 
vectors of the 700-850 hPa layer in each October 
during the 8 year period under consideration 
were analyzed (Figure 8). The 700-850 hPa 

layer is considered as indicative of wind in the lower troposphere, where aerosol transport mainly occurs 
[Dunion and Velden, 2004]. During the second 4 year period (2006-2009), prevailing winds blowing mainly from 
land to sea (Figures 8e-8h) resulted in a drier environment and less precipitation over the east of the Ganges 
basin and northwest BoB (Figure 4a) than during the first 4 year period (2002-2005) (Figures 8a-8d).This caused 
less wet removal of air pollution in the second 4 year period than in the first 4 year period. Second, our analysis 
showed that during the 8 year study period, there was an increasing number of days (Np, in percentage form) in 
each October when prevailing winds blew from land to sea (Figure 9). This suggests some increasing trends in 
the transport of anthropogenic air pollution from their sources in the east of the Ganges basin to northwest BoB. 
Third, for Octobers when Np > 50%, wind convergence was observed over northwest BoB causing the accu- 
mulation of aerosol particles over that region (Figure 10), in line with our previous study [Kishcha et at., 2012]. All 
the three factors contributed to the increasing AOT trend over northwest BoB in the early postmonsoon season. 



Year 

Figure 6. (a-c) Year-to-year variations of AOT for various aerosol species 
(such as sulfates, organic carbon, and desert dust) over zones 6 and 7 in 
northwest BoB in October. The straight lines (dashed for zone 7 and 
solid for zone 6) designate linear fits. 


KISHCHA ET AL. 


©2014. American Geophysical Union. All Rights Reserved. 


1563 


^AGU Journal of Geophysical Research: Atmospheres 


10.1 002/201 3JD020328 



Year 

Figure 7. Year-to-year variations of (a) MERRAero AOT of sea-salt aerosols 
and (b) MERRA surface wind over northwest BoB (zones 6 and 7) in October. 


Figure 1 1 represents maps of the absolute value 
of the magnitude of monthly vertically inte- 
grated mass flux of aerosols including sulfates, 
carbonates, and desert dust over the study 
region. These maps were obtained for each 
October between 2002 and 2009. One can see 
that over northwest BoB (zone 6), the aerosol flux 
was higher during the second 4 year period 
(2006-2009) (when prevailing winds blew mainly 
from land to sea) than during the first 4 year 
period (2002-2005) (when prevailing winds blew 
frequently from sea to land). In October 2009, 
monthly mean wind over zone 6 was minimal 
(less than 1 m/s) (Figure 8 h). As the aerosol flux is 
proportional to wind, then one could expect that 
the aerosol flux over zone 6 should be also min- 
imal. Flowever, as shown in Figure 1 1 h, in 
October 2009, the aerosol flux over zone 6 was 
maximal during the study period. The maximal 
flux under the minimal wind indicates accumu- 
lating aerosol particles over zone 6. This is 
supported by the analysis of wind convergence: 
the wind convergence over zone 6 was also 
maximal in October 2009 (Figure lOe). Therefore, 
we can conclude that vertically integrated aero- 
sol mass flux is sensitive to wind convergence 





Figure 8. (a-h) Spatial distributions of mean MERRA wind vectors of the 700 — 850 hPa layer in each October during the study period. 
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Years 

Figure 9. Numbers of days (in percentage form) in each October during 
the study period when prevailing wind (transporting air pollution) blew 
from the east of the Ganges basin (zone 5) to northwest BoB (zone 6). 


causing aerosol accumulation. Moreover, the 
aerosol flux is sensitive to precipitation: its 
magnitude is maximal over the northwest of the 
Ganges basin and decreases toward northwest 
BoB over the specified zones in the east of the 
Ganges basin, due to wet removal processes 
after significant rainfall (Figure 1 1 ). Therefore, 
the analysis of aerosol fluxes supports our find- 
ings on the effects of wind convergence and 
precipitation on AOT. 

5. Analysis of AOT Increments 

MODIS crosses the study area twice a day. 
Specifically, MODIS-Terra crosses the study area 
at approximately 04:30 UT (10:30 LT), while 








66E 69E 72E 75E 78E 81 E 84E 87E 90E 93E 66E 69E 72E 75E 78E 81 E 84E 87E 90E 93E 


Figure 10. (a-e) Spatial distributions of mean wind convergence (10 6 s ^ of the 700-850 hPa layer for each October, when the number of days with prevailing winds, blowing from the 
east of the Ganges basin (zone 5) to northwest BoB (zone 6), exceeded 50%. MERRA wind reanalysis data were used. 
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Figure 11 . (a-h) Spatial distributions of the absolute value of the magnitude of monthly integrated mass flux (kg m 1 mo ' l ) of aerosols (including sulfates, carbonates, and desert dust) in 
each October during the study period. The absolute value of the flux magnitude was obtained by the following expression: ((Fu) 2 + (Fv) 2 ) 172 , where Fu is the vertically integrated zonal mass 
flux and Fv is the vertically integrated meridional mass flux. 


MODIS-Aqua at approximately 08:30 UT (14:30 LT). MODIS-Terra AOT retrievals were used for updating 
MERRAero AOT at 6 UT, while MODIS-Aqua retrievals were used for updating MERRAero AOT at 9 UT. It is clear 
that after the two consecutive updates, MERRAero AOT at 9 UT corresponds in the best way to available 
MODIS measurements over the study area. During all other times (1 8 h per day), the model simulates air 
pollution over the study area independently of MODIS, using available meteorology and predefined aerosol 
emissions. Taking into account the uncertainty of aerosol emissions over such a complex study area, which 
includes the highly populated and polluted Ganges basin, one could expect some accumulation of model 
errors in AOT simulations during the 18 h period without data assimilation. 

We estimated the NASA GEOS-5 model performance over the study area analyzing model AOT incre- 
ments. These AOT increments are the field differences between MODIS AOT and modeled AOT. AOT 
increments include a complex combination of all model errors in AOT simulations over each specific 
location. To analyze AOT increments, 3 h MERRAero AOT data were used. These 3 h data allowed us to 
distinguish between 3 h periods with and without MODIS AOT assimilation. Figure 12a represents a 
spatial distribution of the 8 year (2002-2009) mean AOT increments at 6 UT. One can see maximal AOT 
increments over the northwest of the Ganges basin (zones 1-3), where precipitation was minimal as 
shown in Figure 4a. In the absence of precipitation, the deficiency in anthropogenic aerosol emissions 
could be the main contributor to the high AOT increments. By contrast, over the east of the Ganges 
basin and northwest BoB (zones 4-7), AOT increments were low, indicating that MODIS AOT assimila- 
tion did not contribute significantly to modeled AOT there. This was despite the fact that in the east of 
the Ganges basin (zones 4 and 5), the population density is the highest in the entire Ganges basin 
(Figure 3), which means the highest industrial and transportation emissions. The most probable reason 
why AOT increments were lower over zones 4 and 5 than over zones 1-3 is the aforementioned 
significant increase in rainfall (accompanied by aerosol wet removal processes) toward northwest BoB. 
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Figure 12. Spatial distributions of the 8 year (2002-2009) mean AOT increments in October at (a) 6 UT and (b) 9 UT. 


Over the east of the Ganges basin and northwest BoB, the NASA GEOS-5 model was capable of 
reproducing the effect of significant rainfall on MERRAero AOT which, in turn, dominates the effect 
of aerosol emissions there. 

Figure 1 2b represents a spatial distribution of the 8 year (2002-2009) mean AOT increments at 9 UT. One can 
see that in spite of the MODIS AOT assimilation at 6 UT, the AOT increments at 9 UT are quite significant over 
the northwest of the Ganges basin (zones 1 and 2). Therefore, MODIS AOT assimilations at 6 UT and 9 UT were 
essential for correcting modeled AOT mainly over the northwest of the Ganges basin, where AOT increments 
were maximal. Over the east of the Ganges basin and northwest BoB where AOT increments were low, MODIS 
AOT assimilation did not contribute significantly to modeled AOT. 

Our findings are further illustrated by Figure 1 3 which represents a comparison between spatial distributions 
of 8 year (2002-2009) AOT trends at 6 UT with and without AOT assimilation. Flere the modeled AOT without 
MODIS assimilation was obtained as the field difference between MODIS AOT and AOT increments at 6 UT. 
One can see some similarity in the two distributions of AOT trends, namely, AOT trends over northwest BoB 
exceed those over the Ganges basin (Figure 13). This is evidence that the increasing AOT trends over 
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Figure 13. Spatial distributions of the 8 year (2002-2009) MERRAero AOT trends (characterized by AOT slopes) in October at 6 UT: (a) with 
MODIS AOT assimilation and (b) without MODIS AOT assimilation. 
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northwest BoB were reproduced by the model not only because MODIS AOT assimilation provided us with an 
opportunity to correct the uncertainty in aerosol emissions but mainly because the model was capable of 
reproducing changes in meteorological factors contributing to the AOT trends. 

6. Conclusions 

The recently developed 8 year (2002-2009) MERRAero assimilated aerosol data set was applied to the study 
of AOT and its trends over the Ganges basin and northwest BoB in the early postmonsoon season. MERRAero 
showed that AOT was lower over the east of the Ganges basin than over the northwest of the Ganges basin: 
this was despite the fact that the east of the Ganges basin should have produced higher anthropogenic 
aerosol emissions due to higher population density, increased industrial output, and transportation. This is 
evidence that higher aerosol emissions do not always correspond to higher AOT over the areas where the 
effects of meteorological factors on AOT dominate those of aerosol emissions. 

Based on the analysis of model AOT increments, we have arrived at the important point of the current 
study which is as follows: the decrease in AOT over the east of the Ganges basin toward northwest BoB 
was reproduced by the model not only because of MODIS AOT assimilation but mainly because of the 
model capability of reproducing meteorological factors contributing to AOT. MODIS AOT assimilation was 
essential over the northwest of the Ganges basin, where AOT increments were maximal, indicating the 
deficiency in anthropogenic aerosol emissions used in the model. Over the east of the Ganges basin and 
northwest BoB, AOT increments were low and AOT assimilation did not contribute significantly to 
modeled AOT there. 

In October, in the absence of aerosol sources in northwest BoB, MERRAero showed increasing AOT 
trends over northwest BoB exceeding those over the east of the Ganges basin. Similar AOT trends were 
obtained using MERRAero AOT with and without MODIS AOT assimilation. Various aerosol components 
showed strong increasing AOT trends over northwest BoB. Therefore, using MERRAero AOT, we obtained 
similar results with respect to AOT trends over northwest BoB, as in our previous study based on MODIS 
data [ Kishcha et al., 2012], 

There were a number of meteorological factors contributing to the increasing AOT trends over northwest BoB: 

1 . an increasing number of days in each October when prevailing winds blew from land to sea, resulting in 
an increase in air pollution over northwest BoB; 

2. during the second 4 year period (2006-2009), prevailing winds blowing mainly from land to sea were 
responsible for a drier environment with less precipitation causing less wet removal of air pollution than 
in the first 4 year period (2002-2005); 

3. in October 2009, the vertically integrated aerosol mass flux over northwest BoB (zone 6) was maximal, 
while monthly mean wind was minimal. This indicates accumulating aerosol particles over northwest 
BoB. The wind convergence over northwest BoB was also maximal in October 2009. Therefore, verti- 
cally integrated aerosol mass flux is sensitive to wind convergence causing aerosol accumulation. 

The MERRAero AOT data set allowed us to determine aerosol species responsible for the AOT maximum 
over the northwest of the Ganges basin, which are both natural aerosols (dust particles) and anthropogenic 
aerosols (carbon aerosols from biomass burning and sulfates). In accordance with MERRAero data, most of 
dust particles fell out due to gravitational settling and wet deposition during dust transport from the 
northwest to the east of the Ganges basin. As a result, over the east of the Ganges basin and northwest BoB 
in the early postmonsoon season, aerosols were dominated by anthropogenic air pollution, such as sulfates 
and carbon aerosols. Our analysis showed that sea-salt aerosols did not contribute to the increasing AOT 
trends over northwest BoB: no sea-salt AOT trend was observed in year-to-year variations during the 
study period. 

We used MERRAero aerosol reanalysis over the Ganges basin and northwest BoB, but this global NASA 
product can be used over other places. Our findings illustrate that while analyzing MERRAero AOT over other 
places, model AOT increments could be helpful in determining (a) areas where effects of aerosol emissions on 
AOT dominate those of meteorology, and MODIS AOT assimilation would be essential for correcting modeled 
AOT; and (b) areas where effects of meteorological factors on AOT dominate those of aerosol emissions, and 
MODIS AOT assimilation would not contribute significantly to modeled AOT. 
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